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1. A Revolution for Lattice QCD? 



2. Improved Actions 



In this lecture I discuss the prospects for sim- 
ulating QCD accurately on very coarse lattices. 
This subject is very important because the cost 
of a Monte Carlo simulation of QCD is very sensi- 
tive to the lattice spacing. Simulation cost varies 
as 
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where a is the lattice spacing, L is the length of 
one side of the lattice, is the pion's mass in 
the simulation, and 5 and e are exponents be- 
tween and 1 that are related to algorithm per- 
formance. With the best algorithms the cost of 
a full QCD simulation QCD increases like 
as a is decreased, making lattice spacing by far 
the most important determinant of the cost. The 
very large exponent means that we that we should 
do everything we can to keep a as large as possi- 
ble. 

Until recently most of us felt that lattice spac- 
ings of order .05-. 1 fm or less were necessary be- 
fore accurate simulations became possible. In this 
lecture I show new Monte Carlo results from sev- 
eral groups indicating that accurate results can 
be obtained even with lattice spacings as large as 
.4 fm. Given that the cost is proportional to 
an a = .4 fm simulation costs lO'^-lO^ times less 
than one with conventional lattice spacings. This 
is a very dramatic improvement — far larger, for 
example, than the difference between a high-end 
PC and a typical supercomputer. If the coarse 
lattices really work, lattice QCD is about to be- 
come vastly more productive. 



A major obstacle to using large lattice spacings 
is that traditional discretizations of QCD are not 
very accurate when a equals, say, .4 fm. We need 
improved discretizations for actions and opera- 
tors. There are two standard approaches to the 
design of improved actions and operators. We ex- 
amine each of them, first for classical field theory 
and then taking account of quantum effects. 

2.1. Perfect Fields (Classical) 

The first approach to discretizing a classical 
field is based upon standard finite-difference ap- 
proximations from numerical analysis; in the con- 
text of lattice QCD this approach is frequently 
associated with Symanzik's workQ. In this ap- 
proach the lattice theory is designed so that the 
lattice field (pi at node Xi equals the corresponck 
ing continuum field (l){x) evaluated at that node:cl 



(3) 



Given this identity the lattice version of nonlinear 
terms in the field equations is exact: for example, 



\(j){x) 



(4) 



Continuum derivatives, however, must be approx- 
imated. Using finite differences, for example, 



A(2) 
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^The analogue of this equation for the QCD gauge field is 

(2) 



Up^ {x) = P exp i I gA ■ dx 



*Invited talk at Lattice '95, Melbourne, Australia (1995). 



where Afj, is the continuum vector potential and P denotes 
path ordering. 
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where 



•H+l 



2(j)i + (f)i. 
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The expansion of d'^(j){xi) is an infinite series in 
powers of a^. Normally only a couple of terms 
are needed. For example, if we are modelling 
a bump in (f>{x) that is 3 or 4 lattice spacings 
across, the term in this expansion is typically 
only of order 10-20% the leading term, while the 
a* term is only a few percent. The terms of or- 
der and higher would be unnecessary for most 
applications in lattice QCD. 

When one wishes to reduce the finite-a errors 
in a simulation, it is usually far more efficient 
to improve the discretization of the derivatives 
(by including more terms in the expansion) than 
to reduce the lattice spacing. For example, with 
just the first term in the approximation to d'^cj), 
cutting the lattice spacing in half would reduce a 
20% error to 5%; but the cost would increase by 
a factor of 2^ = 64 in a simulation where cost goes 
like On the other hand, including the 

correction to the derivative, while working at the 
larger lattice spacing, achieves the same reduction 
in error but at a cost increase of only a factor of 2. 

Although not necessary in practice, it is possi- 
ble to sum all the terms in the expansion of 
a continuum derivative using Fourier transforms: 
we replace 



where 
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'0{\p\<7T/a). (8) 



We do not do this because d^^^ {x-i — Xj ) falls off 



only as l/\xi — Xj\'^ for large 



pf 



resulting 



in highly nonlocal actions and operators that are 
very costly to simulate, particularly when gauge 
fields are involved. So with "perfect fields" it is 
generally far better to truncate the expansion 
than to sum to all orders. However this analy- 
sis suggests a different strategy for discretization 
that we now examine. 



2.2. Perfect Actions (Classical) 

The slow fall-off in the all-orders or perfect 

(2) 

derivative d^f of the last section is caused by the 
abruptness of the lattice cutoff (the 9 funtion in 
Eq. (||)). Wilson noticed that d^^^ can be made 
to vanish faster than any power oi l/\xi — Xj \ if a 
smoother cutoff is introduced |^ . In his approach 
the lattice field at node Xi is not equal to the 
continuum field there, but rather equals the con- 
tinuum field averaged over an interval centered at 
that node: for example, in one dimensional 
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where, in the transform, smearing function 
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provides a smooth cutoff at large p. The lat- 
tice action for the smeared fields is designed 
so that tree-level Green's functions built from 
these fields are exactly equal to the corresponding 
Green's functions in the continuum theory, with 
the smearing function applied at the ends: for 
example. 



ha/2 



ha/2 



dx dy {<jy{x)cj){y)) {12) 

-a/2 Jxj-a/2 
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By rewriting the propagator 



(13) 



in the form 
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we find that the quadratic terms of this lattice 
action are 



5(2) 



\^4>i 4a {Xi - Xj , to) 4>j , 



(15) 



^The smearing discussed here is conceptually simple, but 
better smearings, resulting in improved locality, are possi- 
ble. In [pi the authors use stochastic smearing in which (f>i 
equals tne blocked continuum field only on average. This 
introduces a new parameter that can be tuned to optimize 
the action. 
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where d!^a {xi — Xj , m) is 
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It is easily shown that, remarkably, dpi {x—Xj, m) 
vanishes faster than any power of 1/| 
the separation increases. So the derivative terms 
in such an action are both exact to all orders in 
a and quite local — that is, the classical action 
is "perfect." Again, this is possible because the 
lattice field (j)i is obtained by smearing the contin- 
uum field, which introduces a smooth UV cutoff 
in momentum space rather than the abrupt cutoff 
of the last Section. 

The complication with this approach is that 
it is difhcult to reconstruct the continuum 
fields (j){x) from the smeared lattice fields 4)i. 
Consequently nonlinear interactions, currents, 
and other operators are generally complicated 



functions of the 



and much more difficult to 



design than with the previous ( "perfect field" ) ap- 
proach, where the lattice field is trivially related 
to the continuum field. This is an important issue 
in QCD where phenomenological studies involve 
a wide range of currents and operators. It also 
means that in practice the "perfect action" used 
for an interacting theory is not really perfect since 
the arbitrarily complex functions of that arise 
in real perfect actions must be simplified for sim- 
ulations. (Although recent numerical techniques 
developed by Hasenfratz, Niedermayer and their 
collaborators Q seem to produce accurate approx- 
imations to the classical perfect action even for 
complicated theories like QCD.) 

In general there is no clear choice between the 
two discretization techniques. In the "perfect 
field" approach operators and actions are triv- 
ial to design, but this simplicity may be at the 
expense of locality or accuracy. In the "perfect 
action" approach the action has no a" errors, but 
this perfection requires a very considerable (and 
costly) design effort, even at the classical level. 
As we shall see both approaches have been very 
successful applied to QCD. 



2.3. Quantum Effects 

The techniques described in the previous two 
sections have long been understood, and yet they 
are far from commonplace in lattice QCD analy- 
ses. This is because of quantum effects. Quan- 
tum fluctuations with momenta p> n /a are ex- 
cluded or badly mutilated by the lattice in either 
of the approaches we have discussed. In princi- 
ple we can compensate for this by renormalizing 
the couplings in our quantum opertors. Thus, for 
example, (pd'^cf) is discretized by the replacement 

2 2 

0920^co0A(2)0-ci^0(a(2)) (17) 

where the coefficients cq and ci equal unity in a 
classical theory but are renormalized in a quan- 
tum theory: Ci = I + as^T^ /a)c^^'^ + • • •. The 
rcnormalizations are usually context specific and 
so must be recomputed for each new application. 
However, if a is small enough, p > 7r/a physics 
is perturbative and the rcnormalizations can be 
calculated perturbatively. 

These considerations indicate that improved 
actions are practical only for lattice spacings that 
satisfy the following criteria: 

• a < (important length scales) : This require- 
ment is obvious and applies in classical as 
well as quantum problems. To study the 
static properties of hadrons, for example, 
one cannot use lattice spacings much larger 
than about .4 fm, since light hadrons have 
radii of order .8 fm. On the other hand lat- 
tice spacings much smaller than this should 
not be necessary if improved operators are 
used. 

• perturbation theory works at distances of 
order a or less: The dynamics can be 
perturbative at distances of order the lat- 
tice spacing only if the coupling as{TT/a) is 
small; and if the coupling is small then the 
coefficients in an expansion like Eq. ( [T7[ ) are 
well approximated by their classical values. 
Furthermore the small quantum corrections 
can be systematically computed using (an- 
alytic) perturbation theory. However, if the 
lattice spacing is too large the many ex- 
tra couplings in improved actions and op- 
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erators must be determined nonperturba- 
tively (using, for example, the Monte Carlo 
Renormalization Group), making improve- 
ment very costly. 

The second of these criteria has long de- 
layed the widespread introduction of improve- 
ment techniques into lattice QCD. Until recently 
it was generally believed that perturbation the- 
ory was useless except at very small distances 
(< .1 fm); and therefore improved actions at, 
say, a= .4 fm were also useless, or, at best, pro- 
hibitively costly to generate. This view changed 
completely a few years ago when Paul Mackenzie 
and I noticed that perturbation theory, done cor- 
rectly, works very well even at distances as large 
as 1/2 fmg. 

Perturbation theory is indeed almost useless for 
lattice calculations unless one deals with the "tad- 
pole problem." The nonlinear relation between 
the link field J7^i(a;), from which lattice opera- 
tors are built, and the continuum vector poten- 
tial Afj,{x) leads to tadpole diagrams in lattice 
perturbation theory that do not arise in contin- 
uum calculations. These tadpole contributions 
result in large renormalizations — often as large 
as a factor of two or three — that spoil naive 
perturbation theory. However tadpole contribu- 
tions are generically process independent and so 
it is possible to measure their contribution in one 
quantity and then correct for them in all other 
quantities. Mackenzie and I developed a two-part 
prescription for doing this 1^ : 

• renormalized perturbation theory: A renor- 
malized coupling constant like ayiq*) 
should be used in perturbative calculations, 
not the bare coupling aiat. It is impor- 
tant that this coupling be measured in a 
simulation, for example, by extracting its 
value from that of the expectation value 
of the plaquette; and it is important that 
the scale q* be chosen appropriately for the 
quantity being studied (see ^ for details) D 
By measuring the coupling we automati- 
cally include any large renormalizations of 

*In practice, this procedure is significantly more accurate 
than simply boosting the bare coupling by dividing by the 
plaquette. 



the coupling due to tadpoles. 

• tadpole improvement: Each (x) in a lat- 
tice operator should be replaced by 

U,{X) EE (18) 
Uq 

where uq is the scalar mean value of the 
link, defined to be the fourth root of the 
expectation value of the four-link plaque- 
tte (as measured in the simulation). The 
Mo cancels tadpole contributions, making 
lattice operators far more continuum-like in 
their behavior. 

This prescription for removing tadpole contri- 
butions is now widely used, and its success is well 
documented. I give here only two illustrations. 
In Figure Q I compare (nonperturbative) Monte 
Carlo results for the critical quark mass mda) 
with results from first-order perturbation the- 
ory Q . The critical quark mass is the bare mass 
in Wilson's quark action for which the pion is 
massless. This mass diverges like 1/a in Wilson's 
theory and so should be quite perturbative. How- 
ever, naive perturbation theory, in terms of aiat, 
fails by almost a factor of two even at lattice spac- 
ings as small as 1/20 fm. Reexpressing pertur- 
bation theory in terms of the renormalized cou- 
pling Q!y((7*) leads to excellent results, especially 
considering that this is only a first-order calcula- 
tion. But even these last results are significantly 
improved by first tadpole-improving the quark 
action, and then computing ttIc using first-order 
renormalized perturbation theory. 

In Figure ^ I compare Monte Carlo results for 
— In 1^2,2 at several lattice spacings with predic- 
tions from third-order (renormalized) perturba- 
tion theory 1^; W^2,2 is the expectation value of 
a 2 X 2 Wilson loop. I plot these results ver- 
sus the value of the coupling constant aviq* = 
2.65/a) used in the perturbation theory. This fig- 
ure shows that perturbation theory is still useful 
even when ay ~ 1/2. Similar results are found 
for all sorts of other quantities [|| . Note that 
the 2x2 loops in these simulations range from 
IQ-^ fm to almost 1.2 fm m size. 

These examples and many others show that 
perturbation theory can usually be trusted out to 
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Figure 1. The critical quark mass for Wilson's 
quark action plotted versus the lattice spacing. 
Monte Carlo results are compared with predic- 
tions from perturbation theory using different 
couplings, with and without tadpole- improving 
the quark action. Results are from quenched sim- 
ulations with Wilson's gluon action and /3's rang- 
ing from 5.7 to 6.3; we take a=.l fm at (3 — 6. 
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Figure 2. Minus the logarithm of the 2x2 Wil- 
son loop versus av{q* = 2.65/a). Nonpertur- 
bative Monte Carlo results are shown together 
with predictions from third-order perturbation 
theory. The dotted lines indicate the possible size 
of fourth-order contributions. Results are from 
quenched simulations with Wilson's gluon action 
and /3's ranging from 4.5 to 18. 



distances of order 1/2 fm. I say usually here be- 
cause there are undoubtably situations or quanti- 
ties for which perturbation theory fails to con- 
verge or nonperturbative effects are unusually 
large. However experience indicates that these 
are the exception rather than the rule. Conse- 
quently our criteria for using improved actions 
are satisfied if we choose lattice spacings of or- 
der .4 fm or less; even classically (tree-level) im- 
proved actions, provided they are tadpole im- 
proved, should give accurate results for such lat- 
tices. 

3. Monte Carlo Results 

In the previous sections we have argued that 
improved actions ought to be very effective on 
coarse lattices. Here we review recent tests of this 
idea. The central questions are: Do the correc- 
tions in an improved operator really work? And, 
are the simulations really faster? 

3.1. Heavy Quarks — NRQCD 

One of the most thoroughly studied actions 
that relies upon perturbative improvement is the 
NRQCD action for heavy quarks In NRQCD 
the basic dynamics is given by nonrelativistic 
Schrodinger theory. 
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while relativistic effects and finite-o are included 
by using correction terms: 



SH = -cv 
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The corrections were designed using the "perfect 
field" or Symanzik approach. The couplings are 
computed using perturbation theory, but only af- 
ter the operators are tadpole improved; almost 
all work to date has used just the tree-level val- 
ues for the couplings (c^ = 1). Consequently there 
are only two parameters in the NRQCD action: 
the coupling g and the bare quark mass M. 

Many of the results obtained from NRQCD de- 
pend upon 6H and so provide nontrivial tests of 
the improvement program. For example, spin- 
splittings in the upsilon family of mesons vanish 
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Figure 3. NRQCD predictions for the spin split- 
tings between the IP states of the upsilon. Sim- 
ulations used lattices with a~^«2.4 GeV. 



if is omitted. In Figure 1 1 show NRQCD re- 
sults for the splittings between the IP states of 
the upsilon [|| . The results are in excellent agree- 
ment with experiment, particularly when extrap- 
olated to the correct number of light-quark fla- 
vors, Hf = 3. Note that all of the parameters in 
NRQCD are tuned from the spin-averaged spec- 
trum, and so there are no free parameters in 
this determination of the spin splittings. Note 
also that the spin terms in SH all involve ei- 
ther chromoelectric or chromomagnetic field op- 
erators. These operators are built from products 
of four link operators and so tadpole improvement 
increases their magnitude by almost a factor of 
two at /3 = 6. Without tadpole improvement the 
spin splittings are much too small. 

The correction terms also affect the low- 
momentum dispersion relation of the upsilon. 



Er{p) = Ml 



2 Ah 



(21) 



In a relativistically invariant theory Mi — M2 is 
the upsilon mass. However in a purely nonrela- 
tivistic theory M2 equals the sum of the quark 
masses, which differs from Mi since Mi also in- 
cludes the quarks' binding energy. Only when 
relativistic corrections are added is M2 shifted to 



include the binding energy. Simulations tuned 
to the correct mass Mi = 9.5(1) GeV give M2 = 
8.2(1) GeV when 5H is omitted. When 5H is in- 
cluded, the simulations give Af2 = 9.5(1) GeV, in 
excellent agreement with Mi 

3.2. Gluons — Perfect Field 

The standard Wilson action for gluons has 
finite-a errors of order : 



1 - iRcTrC/pi = Tr(P^,P^,) 



a'^\riRi+r2R2 + nRz 



(22) 



+ 0{a'^) + total derivatives, 
where the are constants and 

R2 - Ep,.,. Tt{D^F,„D^F,„), (23) 
^3 =Ep..,.Tr(i?,,P^,i?,P,.). 

Only operator Pi arises at tree-level; it also 
breaks rotation invariance. The other two op- 
erators are generated by quantum corrections. 
Since operator P3 can be removed by a field re- 
definition. 



(24) 



it doesn't affect spectral quantities and can be 
ignored. 

My collaborators and I have been studying 
"perfect field" (Symanzik) techniques for improv- 
ing the gluon action; most of the new results I 
present in this section are from that collabora- 
tion and are descibed in Following earlier 
work, we added new Wilson loops to the gluon 
action to remove Pi and P2. There are many 
possible choices for the new loops, but among the 
simplest are the rectangle and "parallelogram" : 



C/rt = 



C/pg = 



(25) 



The improved action is |T^,|T^ 

S[U] = /3pi^iRcTr(l-{/pi) 
pi 
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At^iReTr(l-C/,t) 

rt 

/3pg^iReTr(l-C/pg), 



(26) 



with /3pi given as an input, and /3rt and /?pg com- 
puted in tadpole- improved perturbation theory to 
cancel out i?i and i?2- At tree- level, the /3's are 
readily calculated by combining expansions like 
Eq. ( p2| ) for each of the three loops. They are 
tadpole-improved by dividing each Wilson loop 
having L links by (mq)^- One- loop corrections 
have also been computed , but must be ad- 
justed to account for the tadpole improvement. 
We found: 



At = - 



/3p 



'pi 



20 ul 
/3pi 



(1 + 0.4805 a^) 



0.03325 



(27) 
(28) 



Here mq and as are calculated from the measured 
expectation value of the plaquette: 



:iReTr({/p,))^/' 
ln(iReTr([/pi) 
3.06839 



(29) 
(30) 



As in NRQCD, there is no tuning of the couplings 
for the correction terms: tadpole-improved per- 
turbation theory determines them in terms of the 
single bare coupling /3pi. 

The effects of improvement are immediately ap- 
parent if one computes the static quark potential 
using a lattice spacing of .4 fm. Our results for 
both the Wilson theory and the improved theory 
are shown in Figure ^ Off-axis points deviate by 
as much as 40% from a fit to the on-axis points 
for the Wilson theory; these errors are reduced to 
4-5% in the improved theory. The potential from 
the improved theory is essentially unchanged if 
the oig corrections are dropped from the action, 
leaving just the tree-level improved theory (with 
tadpole-improved operators). This suggests that 
0{asa?) corrections, together with 0{a'^) correc- 
tions, are relatively unimportant. Our simula- 
tions also show that the errors are 2-3 times 
larger without tadpole improvement even if as 
corrections are retained. 



Nt 




Tc (MeV) 


2 


6.96 


273 (10) 


3 


7.45 


285 (15) 



Table 1 

The critical temperature for the QCD deconfining 
phase transition computed on lattices with aTc = 
1/Nt using the improved giuon action. The scale 
is set using the IP— IS" splitting in charmonium. 



Morningstar and Peardon, in their presentation 
to this conference, show the first glueball calcu- 
lations with this improved action. They find a 
dramatic reduction in the lattice-spacing depen- 
dence of the glueball mass relative to the Wilson 
theory. Also Bliss, Hornbostel and I have com- 
pleted preliminary measurements of the temper- 
ature Tc of the deconfining phase transition in 
gluonic QCD |l^. Results obtained from lattices 
with a equal to l/2Tc and l/3Tc agree to within 
the statistical errors of a few percent, showing no 
sign of finite-a errors (see Table |l|). 

The spin-averaged spectrum of the "0 family of 
mesons is also independent of the lattice spacing 
for spacings of .4 fm and less[0. This is illus- 
trated by the simulation results for the spectrum 
in Figure ||; the scale for each result is set by the 
IP— 15 splitting. The a-independence (scaling) of 
these results is also evident if they are compared 
with the slope of the static potential V evaluated 
at .6 fm, as shown in Figure ^. The results shown 
for the Wilson theory, however, change by 40% as 
a is decreased, showing large errors. 

As final evidence of the quality of the improved 
actions, I show the IS* and IP radial wavefunc- 
tions for charmonium, computed on lattices with 
two different lattice spacings, using these actions 
(Figure 0). The most striking observation from 
these pictures is that the charge radius of the tp 
is almost exactly equal to one lattice spacing on 
the coarser lattice, and yet the results from the 
coarser lattice are essentially identical to those 
from the finer lattice. These data show that in 
general one needs only a few lattice points per 
direction within a hadron to obtain accurate re- 
sults (few percent) from an a^-accurate action. 
Numerical experiments with the NRQCD action. 
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Figure 4. Static-quark potential computed on G"* lattices with a«0.4 fm using the /? = 4.5 Wilson action 
and the improved action with /3pi = 6.8. 
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Figure 5. S, P, and D states of charmonium 
computed on lattices with: a = 0.40 fm (im- 
proved action, /3pi = 6.8); a = 0.33 fm (improved 
action, /3pi = 7.1); a — 0.24 fm (improved ac- 
tion, Ppi = 7.4); and a — 0.17 fm (Wilson ac- 
tion, /3 = 5.7, from NRQCD Collaboration). The 
dashed lines indicate the true masses. 



Figure 6. IP— 15 charmonimn splitting divided 
by the square root of the slope of the static-quark 
potential evaluated at .6 fm. Results are plotted 
versus lattice spacing for both the Wilson and 
improved gluon actions. 



where corrections are easily turned on and off, 
suggest that 



0(0^) corrections « 15-20% 
0{asa^,a'^) « 2-4% 



(31) 



3.3. Gluons — Perfect Action 

DeGrand, Hasenfratz, Hasenfratz and Nieder- 
meyer have recently computed the classically per- 
fect gluon action for QCD using analytic tech- 
niques for terms quadratic and cubic in and 
numerical techniques for the rest of the action [Q . 
This action is too complicated to use in sim- 
ulations but they have designed a much sim- 
pler, eight-parameter action that approximates 
the perfect action to within a few percent for any 
field configuration. This approximately perfect 
action is built from the plaquette and parallelo- 
gram loops and powers of these. (The powers are 
important.) 

To test their action these authors computed the 
phase-transition temperature Tc for lattice spac- 
ings ranging from l/6Tc to l/2Tc using their ac- 
tion and using the Wilson action. They com- 
pare Tc with the torelon mass a (that is, the string 
tension on lattices of a fixed spatial size). I show 
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Figure 7. The radial wavefunctions for the 15 and IP charmonium computed using improved actions and 
two different lattice spacings. Wavefunctions from a continuum quark model are also shown. Statistical 
errors are negligible for the 15 wavefunction. 
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Nt Wilson Perfect Action 
~ 3.38 (4) 3.09 (3) 

3 3.15 (5) 3.05 (5) 

4 2.98 (3) 3.08 (4) 
Table 2 

Twice the ratio of the square root of the torelon 
mass to the critical temperature for the Wilson 
action and for the approximately perfect QCD 
action. The lattice spacing is 1/NtTc and the 
torelons are for lattices with 2Nt sites in the spa- 
tial direction. 



some of their results for in Table § Re- 

sults for the Wilson action show finite-a effects 
of order 10%; results from their improved action 
are a-independent to within statistical errors of 
a couple percent. There is insufficient data from 
other sources to know whether their results are 
the correct a = results. 

They also computed the static-quark potential 
for two different lattice spacings. The results, in 
Figure |[ show large-r behavior that is indepen- 
dent of the lattice spacing. At smaller radii, re- 
sults from the coarser lattice show some failure 
of rotational invariance. These finite-a errors are 
roughly twice as large as those we find with the 
Symanzik-improved gluon action of the previous 
section. DeGrand and collaborators suggest that 
the difference is because they used the blocked 
gluon operator to compute the potential, which 
introduces additional finite-a errors. I am not cer- 
tain this is the explanation since the errors in 
their measurement are rotationally invariant, and 
the a* errors ought to be pretty small. Another 
possibility is that the errors in the potential are 
due to tadpole effects. The finite-a errors in their 
potential are almost identical to what we obtain 
from the Symanzik-improved action if we omit 
tadpole improvement (with or without as cor- 
rections). This suggests that a tadpole- improved 
version of the perfect action should be tried. 

3.4. Light Quarks 

All standard relativistic quark actions are built 
using a first-order centered difference operator. 



(ignoring gauge fields). Consequently the effec- 
tive lattice spacing in light-quark actions is UeS = 
2a, and not = a as in gluon and heavy-quark 
(NRQCD) actions. This makes the light-quark 
actions the least accurate component in lattice 
simulations and the most in need of repair. Based 
on our experience with NRQCD (Eq. (^), we 
expect order corrections to a light-quark ac- 
tion to be of order 20% when a = .4 fm, which 
corresponds to an UeS of order the light-hardron 
radius. To get errors down to a few percent with 
such a coarse lattice requires an action corrected 
for 0{a^) effects and possibly also for 0{a^). 

The most accurate light-quark action in com- 
mon use today is the tadpole-improved version of 
the Sheikholeslami-Wohlert (SW) action: 

Ssw = ^V'j-A^^) -7-1-771 

where parameter r = 1 is usually chosen. There 
is now considerable evidence showing that this 
action, which has 0{a^) errors, is much more ac- 
curate than the Wilson action for quarks, which 
has 0{a) errors. At this conference Brian Cough 
presented data from Fermilab on quark masses 
that show a significant reduction in finite-a errors 
with the SW action. Also Hugh Shanahan pre- 
sented new UKQCD spectrum results obtained 
using this action. I compare his p mass with 
masses computed by Weingarten and collabora- 
tors [n5l using the Wilson quark action in Fig- 
ure y. I use the QCD scale parameter Ay deter- 
mined from the_measured value of ay (3 A/ a) as a 
reference scale Jj Results from the Wilson quark 
action show large finite-a errors, while the SW re- 
sults are a-independent to within (large) statis- 
tical errors. The two calculations agree within 
errors on the a = mass. Note that tadpole im- 
provement almost doubles the coefficient of cr • F 

^Some people prefer to use the string tension or other 
infrared quantities to fix relative scales. The problem with 
that approach is that most such quantities have significant 
errors themselves, unless they are highly corrected as 
are the NRQCD charmonium and upsilon spectra. Ay, on 
the other hand, when defined in terms of the plaquette has 
no a" errors. The only errors in Ay are from higher-order 
perturbation theory, and these are negligible. 
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Figure 8. The static quark potential computed using an approximately perfect gluon action for two 
different different lattice spacings. 



in the SW action and so is very important to the 
success of this action jl^ . 

These results show that the SW quark action 
is probably sufRciently accurate for most of our 
needs provided a is less than about 1/6 fm (/? = 5.7 
or higher for quenched simulations). However the 
a? errors in this action become noticable at larger 
lattice spacings. This is clear in the SCRI results 
presented at this conference by Robert Edwards 
for the p mass compared with the string tension. 
His results show finite-a errors of order 20-25% 
at a — A fm, which is consistent with my estimate 
inspired by NRQCD analyses. Alford, Klassen 
and I noticed that the leading correction to SW, 
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violates relativistic invariance. We test for this by 
computing the square of the speed of light using 
the pion's dispersion relation: 



Figure 9. The p mass divided by the ay scale 
parameter Ay versus a. Results are shown for the 
Wilson and SW quark actions. The results are for 
quenched simulations with /3's ranging from 5.7 
to 6.2; the lattice spacing is assumed to be .1 fm 
at (3^6. 
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We find that is 30-40% too small when we use 
the SW action with q = .4 fm[p7|. 

The SW action must be replaced by a more ac- 
curate action before useful simulations involving 
light quarks are possible at a = .4 fm. Alford, 
Klassen and I are examining a variety of such ac- 
tions. One of the first we studied is the "D234 
action:" 

<S'd234 = i/i I -A*^^^ ■ 7 -I- m 



-f^A(3).7+^A(4)|^, (36) 

where Ajf' is a discretization of D^, and A^**^ 
is a discretization of ^^D*. This action made 
the error in three or four times smaller, and 
had several other successes as well . I worry 
that quark doublers introduced by the higher- 
order differences may be a problem for this ac- 
tion, but there are several ways of dealing with 
this problem that we are exploring. 

I have no doubt that actions substantially more 
accurate than the SW action will be available 
shortly. The actions discussed above are all of 
the "perfect field" or Symanzik variety. Work has 
also begun on "perfect actions" for the quarks, 
but is still in its infancy ||l8| . 

3.5. Computing Cost 

The use of improved actions on coarse lat- 
tices should lead to substantial reductions in sim- 
ulation costs. My own experience is that the 
reductions are spectacular. Almost all of the 
a — .4 fm results for gluons and heavy quarks 
that I discuss in this paper were generated on 
an IBM RS6000/250 workstation, which performs 
like a mid to high-range PC. These simulations 
typically involved lattices with only 6 points on 
a side, which fit easily in the 8-12 MB of mem- 
ory common on PC's. Quantities like the charmo- 
nium wavefunction can be computed from scratch 
in an hour or two on such a workstation; all previ- 
ous work on charmonium was done with nothing 
smaller than a CRAY, and turnaround times were 
generally quite a bit longer than an hour. 

Beyond the savings associated with smaller 



grids, a large part of the cost reduction in these 
simulations is due to the near absence of criti- 
cal slowing down on coarse grids. Field configu- 
rations decorrelate after 10-20 sweeps even with 
slow algorithms like Metropolis. Note also that 
the asymptotic exponential behavior in propaga- 
tors usually sets in after only one or two lattice 
spacings; this is partly because operators on a 
coarse lattice are automatically smeared (there 
are no high momentum states), and partly be- 
cause the asymptotic region occurs at distances 
of order 1/2-1 fm, which is only one or two lat- 
tice spacings on a coarse lattice. The added cost 
of simulating with a more complicated action is 
negligible (factor of two or three) compared with 
these savings. 

4. Conclusions 

The results I have reviewed here argue persua- 
sively, I think, that classically improved actions 
can be made accurate to a few percent for lat- 
tice spacings as large as .4 fm. Furthermore clas- 
sical (i.e., tree-level in perturbation theory) im- 
provement is sufficient for this purpose, provided 
lattice operators are tadpole-improved. Without 
tadpole improvement, finite-a errors are only half 
removed. 

Some of the features of the gluon actions I have 
discussed were anticipated in actions obtained us- 
ing numerical techniques like the Monte Carlo 
Renormalization Group (see work by Gupta, 
Iwasaki and others). But the high cost of these 
techniques greatly limits their general utility. 
Renormalized perturbation theory and tadpole 
improvement allow us to realize the same goals 
but with negligible costs. 

Both gluon and heavy-quark actions have been 
used with great success at .4 fm, and work is pro- 
ceeding on the development of new, more accu- 
rate actions for the light quarks. Even without 
new light-quark actions, there is now evidence 
that the Sheikholeslami-Wohlert quark action is 
sufficiently accurate for state-of-the-art work at 
a = 1/6 fm (/3 = 5.7 quenched QCD), which is 
already a major advance. 

As we shift to coarse lattices, numerical QCD 
will be redefined by a lO'^-lO^ speedup in sim- 
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ulations. This will have far reaching implica- 
tions for both the sociology of our field and the 
physics it produces. The cntry-lcvcl computer 
needed for serious QCD simulations is no longer 
a CRAY; a personal computer or desktop com- 
puter will do lattice QCD is now accessible to 
almost anyone. The cycle time for ideas is re- 
duced from months to hours — we can dare to be 
much more adventuresome than we have in the 
past. In the short term we should be capable 
of an abundance of 10% accurate results for spec- 
tra, form factors, decay constants, and all sorts of 
other phenomenologically relevant quantities. By 
using our new techniques in simulations on our 
largest computers, we should finally obtain reli- 
able callibrations of the various approximations 
we make in lattice QCD (especially quenching). 
By combining rapid simulations on coarse grids 
with analytic methods, like strong coupling ex- 
pansions, we may move closer to analytic solu- 
tions of low-energy QCD — that is, nonperturba- 
tive QCD without computers. 
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